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Neutrinoless double beta decay {(i(3Qv) is among the only realistic probes of Majorana 
neutrinos. In the standard scenario, dominated by light neutrino exchange, the process am- 
plitude is proportional to m ee , the e — e element of the Majorana mass matrix. Naively, 
qq \ current data allows for vanishing m ee , but this should be protected by an appropriate flavor 

symmetry. All such symmetries lead to mass matrices inconsistent with oscillation phe- 
nomenology. I perform a spurion analysis to break all possible Abelian symmetries that 
guarantee vanishing j3(3()v rates and search for minimally allowed values. I survey 230 bro- 
q , ken structures to yield m ee values and current phenomenological constraints under a variety 

' of scenarios. This analysis also extracts predictions for both neutrino oscillation parameters 

and kinematic quantities. Assuming reasonable tuning levels, I find that m ee > 4 x 10 -6 eV 
at 99% confidence. Bounds below this value might indicate the Dirac neutrino nature or 
the existence of new light (eV-MeV scale) degrees of freedom that can potentially be probed 
elsewhere. 
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I. INTRODUCTION 



Neutrino oscillation experiments have given conclusive evidence that neutrinos have mass and 
' mix. This constitutes the first terrestrial evidence of physics beyond the Standard Model (SM) 

CN ■ and leads to many important questions. For a recent review of neutrino physics see [H, H, El- 

Broadly, it is puzzling why the neutral lepton sector is so different from the other SM fermions 
in both mass scale and mixing pattern. The resolution to this mystery has implications for both 
QO | particle and astrophysics and provides deep theoretical insight into the nature of other high scale 

phenomena. The neutrality of the neutrino under the only unbroken gauge symmetry is the likely 
key to this problem, as it offers the possibility that the neutrino is its own antiparticle via direct 
coupling within a Majorana mass term. Such a mass, as opposed to the more common Dirac mass, 
is composed of only one field and violates all non-zero quantum numbers by two units. The charge 
conjugation properties of the neutrinos, their Dirac vs Majorana nature, are currently unknown 
and their determination is arguably the most important task facing the neutrino community. 

The favored means of probing Majorana neutrinos is via the process of neutrinoless double beta 
decay ((3(30u) where, within a nucleus, two neutrons decay into two protons with no neutrinos 
[3] . This process violates lepton number by two units and may proceed via the virtual exchange 
of Majorana neutrinos. In this case, the decay amplitude is directly proportional to the mass 
of the exchanged electron-type neutrino, or more precisely the e — e element of the Majorana 
neutrino mass matrix (m ee ) taken in the flavor basis where the charged lepton masses are diagonal. 
Current experimental limits on the 76 Ge isotope constrain the f5[5Qu half-life below ~ 10 25 years, 
corresponding to m ee < 0.35 eV at 90% confidence* [1, 0, Eo, U]. Next generation experiments are 
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The translation between measured half-life and m ee is not straightforward, as it depends critically on isotope 
dependent nuclear matrix element calculations, where uncertainties currently range within a factor of three @,@,0|. 
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poised to extend this reach by roughly an order of magnitude to m ee < 0.05 eV [id. Il2l. Il3l|. 

Of course, other exotic interactions can mediate /3/30V, but it was argued in the Blackbox theorem 
that any such Lepton Number Violating (LNV) process will necessarily yield a Majorana 
neutrino mass at some order in perturbation theory, just as a Majorana neutrino mass term will 
lead to LNV processes. This notion was extended in [l6|| to the realistic three neutrino system. The 
authors showed, among other things, that there exists a one-to-one relationship between LNV rates 
and elements of the Majorana neutrino mass matrix such that, in particular, m ee = -4=>- TppQ V = 
and m ee / V pp av / 0. Additionally, using general symmetry arguments, they demonstrated 
that there exists a non-trivial relationship between various mass matrix elements, implying a finite 
set of textures with vanishing m ee . None of these are consistent with the observed oscillation data. 
This leads naturally to the conclusion that if neutrinos are Majorana particles, (3(30u must occur 
at a nonzero rate. 

Here, using similar logic, I explore exactly how small m ee can be under a variety of circumstances. 
If light Majorana neutrino exchange is the only mode of (3(30u, this can be applied directly to the 
interpretation of experimental results. The situation is not as straightforward in the face of other 
contributions, as these will generally effect the mass and LNV rate differently Still, the one- 
to-one correspondence between m ee and TppQ V adds confidence to the conclusion that m ee is a 
good measure of the /3/30f rate for small m ee , even in the presence of arbitrary new physics. This 
statement becomes exact in the limit of vanishingly small m ee . 

This analysis is conducted within the framework of Abelian flavor symmetries acting within a 
three light Majorana neutrino system. It will become clear that the results can be extended beyond 
this paradigm to include non-Abelian groups. Additionally, by the above argument, there is reason 
to believe that the qualitative existence of a lower Tpp 0l/ bound, as well as its connection to m ee , 
should remain valid in the face of arbitrary new physics at scales down to approximately 100 MeV. 
Beyond this point, new light degrees of freedom can contribute directly to the (3(30^ system and 
restrict the result validity. Thus, limits that fall below extracted minimum values are evidence for 
new light physics or the Dirac neutrino nature. A measurement near the derived lower bound would 
indicate a slightly broken symmetry mechanism at work. 

It is natural to wonder how small Tppo,, can be without the introduction of these flavor symmetry 
suppressions. It is well known that, given current neutrino data, there is a well-defined range of 
allowed m ee values. This depends on the neutrino mass hierarchy and oscillation parameters within 
the light Majorana neutrino exchange hypothesis. [HI, Ei, 19, 20, 21|. Neutrinos with normal mass 



spectra can yield vanishing m ee provided appropriate phase and parameter choices. Next generation 
experiments will probe the quasi-degenerate and inverted hierarchy region of the allowed range U 
Clearly, a positive measurement at this relatively large level would not indicate a flavor suppression 
of any kind. The question becomes more involved as smaller values within the normal mass hierarchy 
are explored. When is small m ee no longer an accident? To answer this question, it is instructive 



to take the structure free limit and consider the neutrino mass matrix anarchy hypothesis [22. l23|. 
Here, the underlying neutrino model is sufficiently complicated such that the low energy mass 
matrix appears random and must be treated statistically. In other words, the flavor basis is some 
random rotation from the mass eigenbasis. An analysis of the allowed mixing angle distribution 
is straightforward and described in [2jy]. However, there is an added level of ambiguity introduced 
whenever mass values are discussed, due to freedom in assigning an integration measure to the 
probability distributions. These issues were studied in [22j, l25j|. The distribution of m ee values 



This is likely to improve within the next several years. 



3 



within the anarchy scenario was surveyed in [25l | under a variety of conditions to conclude that 
m ee < 5 x 10 -3 eV implies the existence of a flavor symmetry mechanism, new light degrees of 
freedom, or the Dirac neutrino nature. A (3(30i/ measurement above this limit could be attributed 
to either a flavor symmetry or random fluctuations of the neutrino mass matrix. Below the anarchy 
bound, the present analysis sets a limit on minimum m ee values and identifies what symmetries are 
responsible for the suppression. Above the anarchy bound, it selects those broken symmetries that 
are allowed by the data and makes predictions for other observables that can further constrain the 
system. 

This paper is organized as follows. In Section [III I review Majorana neutrino masses and the 
status of current neutrino data used as constraints in the remainder of the analysis. I then introduce 
flavor symmetries in the context of the neutrino mass matrix and motivate the utility of Abelian 
groups as an ideal laboratory for a comprehensive search for minimal Tppo v rates. In Subsection III Al 
I exhaustively enumerate the flavor symmetry structures that lead to exact m ee = and explore 
their consequences for neutrino oscillation phenomenology. These results are summarized in Table 
HTl which illustrates that of the eleven possible symmetry classes, none are consistent with current 
data. I break these flavor symmetries in Subsection III Bl with the introduction of a single spurious 
U(l)f charged scalar field that acquires a real vacuum expectation value (vev). Subsubsection 
III B II forms the bulk of the analysis, where I numerically survey a comprehensive set of 230 broken 
symmetry structures. For each case (referred to loosely as models), I determine current constraint 
from data, extract the minimally allowed m ee values, and make predictions for future neutrino 
experiments. Variations of these results, subject to improvements in future oscillation parameter 
measurements, are also studied. I conclude in Section II III with a summary of the results and a 
discussion of the limitations of the analysis in the face of new physics. 



II. ABELIAN FLAVOR SYMMETRIES 

If lepton number is violated by physics at some high scale A, the effective low energy Majorana 
neutrino mass Lagrangian term may be written as 

C v = ■^m a pv s ^p- (II. 1) 

In the weak interaction basis where the charged leptons are diagonal, the Greek subscrips are flavor 
indices that run over the three generations e, fi and r. The symmetric mass matrix m aj g may be 
diagonalized to yield positive real mass eigenvalues mi, ni2 and ms by the neutrino mixing matrix 



in the PDG parametrization 261] 




C12C13 S12C13 
U = I -S12C23 - ci 2 S23Si3e i5 C12C23 - si2S23Si3e J<5 S23C13 | e llp2 I (II. 2) 

Sl2«23 - Cl2C23Sl 3 e J<5 -C12S23 - Si 2 C 2 3Si 3 e l5 

that describes the rotation from the flavor basis to the mass basis. I use the shorthand Cjj = cos 9y 
and Sij = sin 6ij for notational convenience. Due to symmetries of the mixing matrix, the mixing 
angles may be constrained within 9ij G (0,7r/2) [271. l28fl with the Majorana and Dirac phases <j>j G 
(0, 7r) and 5 G (0, 2ir) without loss of generality [23,123]. By convention^, the neutrino eigenstates 



^ See for example [28l ] for a summary of mass naming conventions and their relationship with other mixing param- 
eters. 
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Name 


Parameter Combination 


Value 


lcr Uncertainty 


Am^ 


? ? 


1 .UfJ >\ -LVJ C V 


99 x 1 0~ 5 eV 2 




1 777 2 — 777,o 1 
1 7/1,3 '"21 


2.40 x 10~ 3 cV 2 


0.12 x 10~ 3 eV 2 


sin s 


Sin 012 


0.551 


0.017 


sin 9 a 


sin 023 


0.707 


0.046 


sin R 


sin 13 


0.1 


< 0.14 









0.50 eV 


£ 


toi + m 2 + m 3 





0.24 eV 




\micl 2 cf 3 + ?n 2 s?2 c i3 e2i02 + m 3 .s2 3 e 2j ^ 3 - d )| 





< 0.175 



TABLE I: Summary table of current neutrino results together with naming conventions and parameter 
definitions. Columns three and four list best fit central parameter values and lcr uncertainties. The first 
five entries were adapted from the global oscillation analysis of jH3|. The upper bounds for the last three 
entries come from endpoint measurements of Tritium /3-decay (wi„ c ), cosmological observations (£) and 
(3(3Qv, respectively [llj. These are used as model constraints through the remainder of this analysis. 



are ordered in mass squared such that m 2 , and m 2 . have the smallest separation, the so-called solar 
mass squared splitting Am 2 s = m\ — m\ > 0, while m.3 is the most distant state. The identity 
of the lightest eigenstate depends on the mass ordering, with v\ and 1/3 existing as the lightest 
state for the normal and inverted hierarchies, respectively. Neutrino oscillation data constrains the 
mass eigenvalue squared differences and mixing angles as can be seen from the first five entries 
of Table [p which lists both the current best fit values and la uncertainties, as adapted from [30]. 
See also [ll|, [3l|, |32j]. Additionally, kinematic probes of the endpoint of the tritium beta decay 
spectrum, cosmological observations and even (3(301/ constrain the absolute neutrino mass values 



[101 ] as shown in entries six, seven and eight. To date, these have yet to observe positive signals 
but have been successful in bounding neutrino masses below the eV level Next generation 

experiments will extend the reach of tritium decay and cosmological measurements to m Ue < 0.2 eV 



331 ] and E < 0.1 eV [34j, |3J3, [3fj] , at 90% confidence, respectively. The single Dirac phase 5 and two 
Majorana phases (^2, and ^3 are currently unconstrained by experiment. The first column of Table 
[J lists the parameter name conventions incorporated in this analysis. The subscripts S, A and R 
attached to the oscillation parameters refer respectively to "solar", "atmospheric" and "reactor," after 
the primary /historical neutrino sources used in their measurement. These parameter constraints 
must be satisfied by all viable neutrino mass models. 

The entries of the symmetric matrix m a p are constrained to be small by a combination of the 
kinematic and oscillation neutrino data. The suppression of this term with respect to the other 
charged fermions is likely due to the high scale of new physics. The mass term of Eq. (|H.1|) is the 
relic of an effective operator after electroweak symmetry breaking and as such, m a p oc v (j)" a a p 
where v = 0.246 TeV is the Higgs vev and a a p is a matrix of complex constants. This may arise 
from a simple dimension five operator as in the seesaw mechanisms 37, 38, 39, 40, 41, \^ or from 
more exotic high dimensional interactions [13, |43|, [4J]. The specific UV completion is irrelevant 
for the purposes of this paper. All that is required is a possibly broken flavor symmetry principle 
that is manifest in the low energy effective system. Then the neutrino fields of Eq. (IH.lj) will 
transform under a representation of the symmetry group and only those mass matrices m a p that 
render L v invariant will be allowed. Some of these symmetries will require particular texture zeros 
for full invariance that when broken will induce small deviations from zero that may be probed 
experimentally. 

Assuming the Majorana nature of neutrinos, I extract the minimal value of m ee for a variety of 
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model classes. Abelian flavor symmetries are useful in this endeavor due to the freedom of charge 
assignments to the neutrino fields. Under non-Abelian symmetries, on the other hand, the fields 
will transform under some representation of the group and the Lagrangian terms will be composed 
of invariant field combinations. Members of the field multiplets are assigned quantum numbers, as 
in the Abelian case, except that these are imposed by the group representation. These charges are 
an important factor in mass matrix construction and may be mimicked by a properly constructed 
Abelian symmetry. There is more freedom in mass matrix construction associated with Abelian 
symmetries, which proves to be useful when scanning for small m ee values. The restrictions imposed 
by non-Abelian symmetries can only push the extracted m ee values up. Thus, in what follows, the 
discussion is restricted to Abelian flavor symmetries. 



A. Unbroken Abelian Flavor Symmetries 

Here I assume that a zero m ee element is protected by an Abelian flavor symmetry. Given the 
U(l)f charges n e , and n T supplied respectively to the u e , and v T neutrinos^, it is simple to 
derive the form of the resulting symmetric Major ana neutrino mass matrix invariant under U(l)f 

m a/3 = Ma a/3 5 (na+n ^ !0 . (II. 3) 

a a p is some, presumably 0(1), complex constant and M ~ v (x)™ 1S the neutrino mass scale. 
Currently, M is bounded at 0.05 eV from below by the atmospheric mass squared difference [111. \3l\ 



and from above at roughly 1 eV by cosmological data [U|. Notice that the electron neutrino mass 
term m ee v%v e has charge 2n e . To preserve the imposed flavor symmetry, either n e = or m ee = 0. 
Thus, a non-trivial transformation of v e under U(l)f guarantees m ee = and consequently a 
vanishing f3f30v rate. If one assumes that all allowed entries are nonzero, there are eleven possible 
m ee = mass matrix classes that can be obtained from Eq. (|II.3|) by scanning charge assignments^. 
These are listed in Table IH using the charge assignment notation (n e ,n^,n T ), along with their 
associated mass matrix and neutrino mixing predictions. There are only a small number, between 
one and three, of free parameters in each matrix entry. These must be used to construct three mixing 
angles and two mass squared difference predictions upon diagonalization. Thus, one obtains large 
correlations among the derived oscillation parameters. If the coupling coefficients a a p are allowed 
to take on any value, the class assignments are superfluous in that some entries are just special 
limiting cases of other classes. For example, Cll is just a special case of CI with a w = a TT = 0. 
These distinctions are made here due to the expectation that all matrix elements allowed by U(l)f 
should be of the same order, in which case each class yields different predictions. The predictions 
are obtained by a simple diagonalization of the resulting mass matrix under the convention that 
the smallest mass squared difference defines the solar oscillation frequency and the next largest 
the atmospheric. The largest mass squared difference is the sum of the smaller two and converges 
to Am\ when Am 2 s /Am\ is small, as required by data. Degenerate eigenvalues are treated as if 
they possessed small splittings induced by symmetry breaking effects in anticipation of Subsection 



* The neutrinos are components of the SM left-handed doublet fields L a = ( Ua ) , thus left-handed charged leptons 



la will also be charged under U(l)f. The symmetry dictates the identity of the flavor basis by well defined charge 
assignments. For my purposes, freedom in the flavor charge structure of the right-handed leptons en can be used 
to propertly construct the diagonal charged lepton mass matrix. 

The trivially zero mass matrix, obtained when n a + np for all flavors a and /3, is not included in this listing. 
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III PI The split levels are then associated with the solar mixing sector and interpreted as such to 
make mass hierarchy predictions. Even then, the neutrino mass ordering can only be predicted 
when the lightest eigenvalue vanishes. This occurs in all classes except CI, C3 and C6. In these 
cases, one may derive relationships between the mixing parameters and discrete hierarchy choices. 
Degeneracies leading to invariant matrix subspaces in classes C2 and C4 yield additional freedom 
corresponding to an arbitrary rotation within the invariant subbasis. This is parameterized by 
a mixing angle 9 in Table [U that may take on any value. It should be noted that symmetry 
breaking effects of U(l)f, which select a definite mass basis, destroy this freedom by selecting a 
particular value of 9. Here, the symmetry breaking mechanism yields discrete variable changes and 
is therefore more important than in the other cases where deviations from the predictions of Table 
Ull are parametrically small, or proportional to the symmetry breaking order parameter. The goal 
of this exercise is to understand how close each £7(1) symmetry structure comes to reproducing the 
current neutrino oscillation data, and in what ways they tend to fail. 

Most predictions of Table HT1 are well defined and need little explanation. CI, however, is more 
involved and deserves a separate discussion for clarity. Changing the parametrization of the mass 
matrix for convenience to 



(II.4) 



I find that the lightest mass eigenvalue is zero, while the absolute mass squared differences are 
(a — d) 2 and Aad. The associated atmospheric and solar parameters depend on their relative 
splitting sizes. We are left with two cases defined by their predicted mass hierarchy. For the 
inverted hierarchy, Am| = (a — d) 2 and Am^ = Aad, while all mixing angles vanish except 9r, 
which is given by s 2 R = (d — b)/2d. For the normal hierarchy, the mass squared differences are 
reversed and all mixing angles vanish except 9a 1 which is given by t\ = (d — b)/(d + b). It is 
interesting that this is the only case, due to availability of three free mass matrix parameters, that 
allows for a nonzero solar mass squared difference. Consequently, this class is well-suited to fit 
the neutrino data with only minor modifications. While neither case is consistant with neutrino 
oscillation phenomenology, it is clear that the normal hierarchy choice can be pushed "closer" to 
the observed form. For nearly maximal atmospheric mixing, b must be small. For the mass squared 
differences to work out, the parameters a and d must be unnaturally tuned to a /Ad 10 ±2 . Hence, 
this scenario is far from ideal when considered with universally 0(1) parameter values. This mass 
matrix texture is theoretically motivated by variants of /i — r symmetries and is commonly found 




in the literature. See for example [13, |45|, |46(] and references therein. 

Each entry of Table [III defines a class of models with similar characteristic predictions. It is 
important to note that none of these classes fit the neutrino mixing data. This is another reiteration 
of the fact that for Majorana neutrinos, m ee / 0, which implies a nonzero f3f30v rate. In particular, 
none of the classes predict realistic mixing angles. This can be seen by inspection, as most cases 
predict some combination of maximal 9r or vanishing solar/atmospheric angles. Classes C7 and 
C8 are less trivial but still ultimately fail. C8 fails due to the relationship between the angles 
sr = tstA, implying that a small Or must be accompanied by either a small 0s or 9 a- Similarly, 
in C7 0s must be close to it/ 2 to insure a small reactor angle. Additionally, all classes except the 
first yield degenerate eigenvalues implying Am| = 0, which further contradicts observation. 
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Class 


Charge 


Matrix 


Predictions 


CI 


K,o,o) 


M 


/ \ 

a w a MT 
\ a pT a TT / 




A J. r\ /"\T 1 TT* 1 \ 

sr : = 0, is = (Normal Hierarchy) 
i A = 0, is = (Inverted Hierarchy) 
Can tune A 2 m's to fit data. See text for details. 


C2 


(n e ,0,n r ) 


M 


/ \ 

a MM 
V / 




A 2 n A 2 71 f 2 2 

Am s = 0, Ami = M Sm 
s fi = sin#, t A = cos 9, t s = 
Normal Hierarchy 


C3 


(n e ,0, -n e ) 


M I 

\ 


a eT 

a er 


) 

/ 


A 2 nIA 21 7i/f2l2 21 

Am s = 0, |Ami| = M z \a z TT - a 1 ^ 
4 = 1/2, t A = 0, i s = 
Any Hierarchy 


C4 


(n e , 71^,0) 


M 


/O \ 

nn ° 
\ O a TT j 




A 2 n A 2 71 ^2 2 

Am z s = 0, Am^ = M*a TT 
sr = 0, t A = 0, is = tan 9 
Normal Hierarchy 


C5 


(n e ,n e , — n e ) 


M ^ 


a er 
a^r 


) 


A 2 n A 2 71 r2 / 2 i 2 \ 

Am s = 0, Am A = M z (a z eT + a^ T ) 
s% = 1/2, t A = 0, is = a^r/aer 
Inverted Hierarchy 


C6 


(n e ,-n e ,0) 


M | 


a el j, 

a eM 
a TT 


) 


Am| = 0,|Ami| = M 2 | a 2 r „ a 2j 
s fl = 0, t A = 0, ig = 1 
Any Hierarchy 


C7 


(n e ,-n e ,n e ) 


ill ^ 


a elJ 
a Mr 


) 


Ar4 = 0, &m\ = M\al T + a%) 

2 2 2,2 
„2 %i + 2 +2 a e^+ a M T 

Inverted Hierarchy, s# = t A /ts 


C8 


(n e , -n e , — n e ) 


„( 


Q>efi Q-er 

a e ^i 


) 


Am| = °» Am A = A^Oer + Oep) 
„2 _ a\ T f 2 _ «er +2 _ a L 

Inverted Hierarchy, = i^is 


C9 


(n e , ~n ei n T ) 


M 


/ a e „ 0\ 

a e/J 
\ / 




A 2 n A 2 71 r2 2 

Am s = 0, Am^ = MX,, 
sr = 0,t A = 0,t s = l 
Inverted Hierarchy 


CIO 


(n e ,n. M ,-n e ) 


M 


/ a eT \ 


/ 




A 2 n A 2 7i //"2 2 

Am s = 0, Am^ = M a eT 
4 = 1/2, t A = 0, is = 
Inverted Hierarchy 


Cll 




M 


/O \ 

a MT 
\0 a MT / 




Am| = 0, Am^ = M 2 a 2 r 
4 = 1/2, i A = 0, is = 
Inverted Hierarchy 



TABLE II: Exhaustive summary of neutrino mass matrices with m ee = protected by Abelian flavor 
symmetries. Column one defines the symmetry class as referred to throughout the text. Columns two and 
three list representative U(l)f charge assignments to the neutrino flavor eigenstates and the implied mass 
matrix, respectively. The final column shows oscillation parameter predictions for each class. None of these 
are consistent with current neutrino data. See the text for more information. 



B. Broken Abelian Flavor Symmetries 

Exact m ee = may be excluded by this line of reasoning, but it still may be unobservably small. 
In this case, the symmetries of Table [III should still be approximately valid, only broken by some 
small amount e. These broken scenarios should still retain some of the features of their parent 
classes, such as mass hierarchy or large/small mixing angles. Hence, one would expect that classes 
such as CI, C2, and C4 with multiple predictions that can be made consistant with data, will be 
broken far less than, say, classes C10 and Cll that are far from data. 

I parameterize this symmetry breaking with the introduction of a spurious scalar field s, charged 
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under U(l)f, that acquires a nonzero vev. For this spurion analysis, s is just a mathematical 
construct used to understand the pattern and size of symmetry breaking. Specifically, I assume 
that s has U(l)f charge n s and acquires a real vev e. With this, small mass term corrections are 
induced to help fit the data. The form of the resulting mass matrix may now be described by a 
4-tuple (n e ,n^,n T ;n s ). For example, the symmetry of class C6 may be broken to yield 



/ 2 


a eT e N 




m(n e , -n e , 0; n e ) = M 1 a e)1 a^e 2 


v e 


| , (H.5) 


\ a eT e a^e 












m(n e , —n e , 0; 2n e ) = M 1 a efl a m e 


"') 


(II.6) 


\ 


a TT / 





depending on the spurion charge assignment. In what follows, I will refer to these broken structures 
by integer valued 4-tuples. In the previous example, this would be (1,— 1,0; 1) and (1,-1,0; 2), 
respectively. These are not unique but adequately represent a broken class. In parent classes CI, 
C3, and C5 — C8 with a single assigned charge, the non-trivial^ spurion charges are always multiples 
of n e . These cases yield between two and three perturbed mass matrix structures for each class. 
For the remainder of the parent classes, with two distinct charge values, non-trivial spurion charges 
must be multiples of either neutrino charge value or their sum. This freedom leads to a proliferation 
of perturbed mass matrices - between 37 and 52 for each class. For a general charge assignment, 
the lowest order mass matrix elements are given by 

oo 

m al3 = MY,a%e l 5 {na+np)tma , (II.7) 
i=0 

where a^l are complex constants for the 0(e l ) terms. Once again, it is reasonable to assume that 

(i) 

all nonzero a a n are of the same order of magnitude, since they all arise as coupling constants within 
the spurion included invariant Lagrangian. Clearly, the delta function vanishes for all but a single 
i value. From here, it is trivial to build up the matrix to higher order. If the operator 0% is some 
U(l)f invariant function of the neutrino and spurion fields, of order 0{e % ) after flavor symmetry 
breaking, then the operator = ssOi is also invariant of order 0(e l+2 ) after s acquires a vev. 
This process may be continued indefinitely to yield the full matrix. Thus, each nonzero matrix 
element is expressable as a power series in e 2 . For most cases of interest, where e is sufficiently 
small, all but the leading terms described by Eq. (|II.7|) may be neglected. This new broken matrix 
may be diagonalized perturbatively to yield predictions for oscillation and kinematic parameters as 
a function of e. Due to the number of charge assignments and breaking patterns, it is most efficient 
to study this numerically. 

1. Numerical Results 

I now numerically survey all mass matrices generated from a broken U(l) flavor symmetry to 
order e 10 in m ee and e 3 in all other elements. The unequal treatment of the mass matrix is due to the 



Nontrivial in this context refers to an assignment that yields perturbed mass matrices distinct from the parent 
class structure. 
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search goal of small m ee . Corrections higher than e 3 make little difference to the fit, but one must 
be able to distinguish cases with different large m ee suppressions, even when the rest of the matrix 
is identical. It will turn out that no structure with m ee suppressed by more then e 8 is consistant 
with data. Thus, truncating the search at e 10 is safe. Using the notation of Eq. (|II.7p . I set a ee = 1 
in order to remove the ambiguity of simultaneous rescalings of all the a a p and M. Furthermore, 
I rescale e' — > e gcd ( p \ where P stands for the set of all nonzero powers of e in the original mass 
matrix and gcd is the greatest common denominator function. This removes another unphysical 
ambiguity. To see this, suppose that two different charge assignments yield the same perturbed 
matrix structure up to an overall rescaling of e — > e> . In this case, the best fit to the data would 
select out preferred e values related by the rescaling, but the same mass matrix elements. Thus, 
from the neutrino mass generation standpoint, both charge assignments yield the same predictions 
and one should be selected to represent the system**. After all rescalings and truncations, there 
are 230 distinct possibilities that may be categorized into one of the eleven classes of Table HH In 
each case, I scan the symmetry breaking parameter e, the mass scale M and coupling constants 
a^g to fit current neutrino constraints listed in Table [J and extract the smallest allowed m ee . 
Specifically, for each model I perform a x 2 fit to minimize the function 

e (A2ro; : A2mj>2 + e (3m " ; r m ^ + ^ + gy. (h.b, 

j={S,A} (J A^m j j={S,A,R} a smd J U m v 

where the primed quantities are evaluated from the broken mass matrix and the unprimed best fit 
values and uncertainties are taken from Table [B Here, only M and a£l are varied to yield the one 
parameter function F(e). When F(e) falls below a critical value, the broken flavor symmetry is 
allowed at a specified confidence. To be conservative, I use 99% confidence limits throughout this 
analysis. The critical \ 2 value depends on the number of degrees of freedom within the system 
and is therefore different for each case. The allowed domain where F(e) falls below its critical 
value may then easily be scanned for the smallest m ee value. If no such region exists, the charge 
assignment and breaking structure is disfavored by current data at 99% confidence. To maintain 
the perturbativity of the system, I hold e < 2/3; in which case, the largest possible (fourth order) 
corrections are only (2/3) 4 « 20%. Typically, corrections will be much smaller than this since m ee 
goes like some power of e and small m ee values favor small e. For this procedure to make sense, 
the constants \a a p\ must be constrained by some naturalness criterion, else the small e values may 
be compensated by large coupling constants resulting in the loss of algebraic structure information. 
In this spirit, I only allow the constants to vary symmetrically about unity, in a log 10 sense, by 
a small amount. That is, 10 _c < \a a p\ < 10 c for some small number c. A one and two order of 
magnitude spread is defined by c = 0.5 and c = 1, respectively. Ideally, the relative size of each 
matrix element should be determined by e alone, so it is clear that c should not be much greater 
than 0.5 for a typical e ~ 10 _1 . Even this range is dangerous near the upper e limit where order e n 
terms can easily be larger than order e n_1 terms. To remove this problem completely for the full e 
range, one needs c < 0.1 which offers very little parameter freedom and is not realistic. I present 
data for c = 0.1, 0.3, 0.5, and 0.7 in the attempt to cover a comprehensive range of naturalness 
criterion. 

Figure [T] histograms the 99% allowed flavor structures by minimum m ee value for c = 0.1, 0.3, 0.5 
and 0.7 on a log 10 scale. These are color-coded by parent class. Reference lines indicating the 0/3Qv 



** This is not true when a spurion analysis beyond neutrino mass is performed. In that case, each assignment would 
yield distinct predictions for other processes. 
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FIG. 1: Histogram of minimal m ee values extracted from data allowed models, color coded by model class. 
From top to bottom, the panels represent the c = 0.1, c = 0.3, c = 0.5 and c = 0.7 cases, c defines fine 
tuning sensitivity by limiting the relative magnitudes of the free matrix parameters as 10~ c < \a a p\ < 10 c . 
A larger c value represents more parameter freedom, which leads to a greater number of allowed models and 
smaller m ee , as is clear from the plot. See the text for more details. 



current bound (solid black) and future reach (broken black) are included for reference as well as the 
anarchy bound (solid gray). Due to increasing parameter freedom, the number of allowed models 
increases and broadens as the c value is pushed higher. For example, the 57 models of c = 0.1 
span less than four orders of magnitude, while the 74 models of c = 0.7 span seven. The smallest 
m ee values for c = 0.1, 0.3, 0.5 and 0.7 are 6.43 x 10~ 5 eV, 2.96 x 10~ 5 eV, 3.92 x 10~ 6 eV and 
5.28 x 10~ 7 eV, respectively. Optimally minimized charge assignments along with their relevant 
parameters are summarized in Table IHIl The mean e values for each c is also shown for easy 
comparison. The broken symmetry that yields the smallest rates are (4, —1,0; 1) for the c = 0.1, 
c = 0.3 and c = 0.5 cases and (4,0, —1; 1) for the c = 0.7 case. The last symmetry, optimal for 
c = 0.7, only beats out (4, —1,0; 1) by 0.73%. These are members of very similar parent classes C4 
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c Case 


Class Charge fn ee (eV) e < e > sin. Or sin 9a sin $5 Hierarchy 


0.1 


C4 (4, -1, 0; 1) 6.43 x 10~ 5 0.402 0.398 0.0870 0.5648 0.5341 N 


0.3 


C74 (4, -1, 0; 1) 2.96 x lO" 5 0.376 0.367 0.0863 0.5793 0.5293 N 


0.5 


C4 (4,-l,0;l) 3.92 x 10~ 6 0.276 0.279 0.0859 0.5945 0.5217 N 


0.7 


C2 (4,0,-1; 1) 5.38 x 10~ 7 0.203 0.218 0.0859 0.5932 0.5207 N 
C4 (4,-l,0;l) 5.32 x 10" 7 0.203 0.0862 0.5966 0.5200 N 



TABLE III: Summary of models with the smallest m ee values for the c = 0.1, c = 0.3, c = 0.5 and c = 0.7 
cases, as listed in column one. Columns two, three and four list the model's parent class, a representative 
charge assignment and derived m ee . Column five shows the optimal amount of U(l)f symmetry breaking 
via the order parameter e to be compared with the mean value < e > in column six. The last four columns 
give mixing angle and mass hierarchy predictions. Two items with similar m ee are listed for case c = 0.7 for 
easy comparison. See the text for details. 




and C2. It is not a coincidence that these both contain invariant subspaces that allow for additional 
parameter freedom and predict the normal mass hierarchy. m ee is so small in these cases because 
it is suppressed by e 8 . The reason there is so great a difference between them is due solely to the 
increased parameter freedom of higher c values allowing for smaller e. Even these tiny variations in 
e are amplified in the m ee relations and can easily yield order of magnitude differences. 

Inspection of the class descriptions of Table [TT] reveals that both classes C2 and C4 have two 
distinct problems that must be solved by symmetry breaking. Specifically, it must induce a nonzero 
Am|, as well as push 9$ and 9a, respectively, up to allowed levels. Additionally, the symmetry 
breaking mechanism must be able to explain the 9 values required by data, namely 9 ~ and 9^9$ 
for classes C2 and C4, respectively. Looking at the broken C4 case (4, —1,0; 1) as a representative 
example, it is easy to see how these are solved. An allowed real valued matrix corresponding to the 
c = 0.3 case is 

-1.42e 3 

2.00e 2 0.78e | + 0(e 4 ), 
0.78e 0.50 

where M = 0.074 eV and e = 0.37. The m ee oc e 8 is included in the matrix for illustrative purposes. 
Some predictions of this structure are sin 9 r oc e 3 , sin#s oc e and Am|/M 2 oc e 4 up to 0(1) 
modifications. Taking the allowed Am| and neutrino mass scale range in the last relation implies 
e between 0.01 and 0.25 as observed in the fit. The large solar mixing angle data selects the upper 
part of this range which leads to a large 9r prediction. Once the numerical factors are accounted for, 
this matrix structure is very similar to that of class CI, with the lower right (/, — r block elements of 
the same order, disagreeing at most by 56%. As previously discussed, CI is a popular texture that 
is allowed, provided small symmetry breaking. The charge assignment (1, 0, 0, 1) of class CI fits the 
data in the c = 0.5 scenario with e = 0.0158, yielding the e 2 suppressed m ee = 1.97 x 10~ 5 eV. It 
turns out that the majority of the minimized matrices with m ee < 10~ 4 eV acquire this form. This 
is optimal from the small m ee perspective, since all data can be accommodated with a vanishingly 
small m ee element. That is, the e — e matrix element does not significantly contribute to the fit. 
I find that no broken symmetry with an exact e — e texture zero can yield this approximate CI 
structure, but highly suppressed terms are possible. For these cases, the size of m ee depends most 
critically on the e power suppression as opposed to the amount of symmetry breaking. 

A handful of these, particularly those of class C9, will be probed by next generation (3(30is 
experiments. For example, in the representative c = 0.5 case, the class C9 charge assignment 
(1,-1,-3,1) will be explored. The breaking is optimized with e = 0.38, leading to an e 2 sup- 
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FIG. 2: Scatter plot projections of neutrino oscillation parameter c = 0.5 predictions taken from data 
allowed models. 90% and 99% parameter confidence limits are indicated by the dark and light gray regions, 
respectively. 

pressed m ee = 0.070 eV with the inverted mass hierarchy. It should be noted that null result 
bounds set by (3(3Qv experiments do constrain those models with higher minimal m ee values. How- 
ever, a measurement of (3(30u does not pick out a particular broken flavor model since all that is 
being plotted is the minimum m ee value. In other words, all unconstrained models should be con- 
sidered equal candidates. Furthermore, bounds or measurements of small m ee only indicate a flavor 
symmetry when they fall below the anarchy bound near 5 x 10~ 3 eV. Above that value suppressed 
could simply be a random fluctuation of an anarchical mass matrix. 
Each of these models, taken at their minimum m ee values, make predictions for other observ- 
ables. These are found by diagonalizing the perturbed mass matrix, together with the minimization 
parameters, to obtain the mass eigenvalues and mixing angles. They may also be combined, as pre- 
scribed in Table [U to yield the mass squared differences as well as the effective neutrino mass 
relevant to tritium beta decay and the cosmological neutrino mass sum. This prescription also se- 
lects the neutrino mass ordering predicted by each model. Figure [2] displays scatter plot predictions 
of the c = 0.5 models allowed by data shown as projections onto the neutrino oscillation parameter 
space. The dark and light gray shaded rectangles illustrate the 90% and 99% allowed regions, 
respectively. These are not true confidence contours, as they do not take correlations into account, 
but adequately reflect parameter regions allowed by data for the purposes of this analysis. Most of 
the scattered points are found to possess the normal mass ordering. There are 64 cases predicting 
normal and 7 cases predicting inverted mass orderings. Generally, the normal hierarchy cases yield 
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FIG. 3: Scatter plot projections of neutrino kinematic and absolute mass predictions (m„ c , £ and mee) 
taken from data allowed models. Current 99% confidence limits are indicated by the gray region, and the 
solid black contour line illustrates the next generation experimental reach. 

a lower m ee . A quick inspection of this plot reveals that the mass squared differences, particularly 
Am|, have little spread and are thus not suited to distinguish between the symmetry structures. 
Am\ may be tuned to fit data in even the unbroken cases enumerated in Table [III with minimal 
effects on remaining parameters. Additionally, Am| is a relatively small splitting that arises when 
U(l)f is broken in all classes but CI, where it is present from the beginning. For models that 
allow such small splittings, it is easy to adjust parameters to fit the data. Mixing angles, however, 
turn out to be a much better diagnostic tool. The unbroken classes of Table [TT] predict vanishing 
or maximal angles for the majority of entries. If these fit the data as given, say, maximal 6a and 
vanishing Or, symmetry breaking is likely to drive these away from the preferred values. In the 
opposite limit, when predictions are far from data, the symmetry breaking must drive the angle 
significantly to reach the allowed range. This results in a large spread in predicted angles. I find 
that Or and 9a are the best laboratories for constraining the model classes. The points in the 
projection onto this parameter plane spans almost the entire allowed region. The triangular shape 
of the scatter profile in this panel is an artifact of the fitting procedure. If a large deviation from 
best fit is found in one parameter, there is little room left for deviations in any other parameter. 

In a similar way, Figure [3] displays c = 0.5 prediction scatter plots in the space of absolute mass 
observables m Ue , £ and m ee . The m ee direction is histogrammed in the third panel of Figure [1] 
and is included here to help visualize the relationship among the parameters. Once again, light 
gray illustrates the 99% allowed region. The solid black contour line represents the reach of next 
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generation experiments. These contours only approximate the present and future experimental 
bounds, as there are large correlations among m ee , E and m u under the three light Majorana 



neutrino assumption employed here. See for example [2l|, [47|, |48|, |49j, [5CJ and references therein 
for a discussion of these correlations. The smallest inverted hierarchy m ee = .0223 eV is for the 
C7 charge structure (1,— 1,1; 1). This is the only class CI model consistent with data, and thus 
requires a relatively large amount of symmetry breaking, with e = 0.3159. The remaining six 
inverted hierarchy cases are of classes C9, CIO and Cll, consistent with expectations from the 
Table [III predictions. Within these, the hierarchy determination depends heavily on the specific 
charge structure. Those predicting normal spectra generally have larger symmetry breaking in 
order to overcome the influence of their parent classes. Nevertheless, these can have m ee as small as 
10~ 5 eV due to large e power suppressions. Still, it is gratifying that the inverted cases are all above 
the 99% confidence lower limit defined in terms of the m ee entry of Table [J with 771,3 = and 4>2 = vr. 
Evaluated, this is m ™ ted » cos 2 9 13 cos A ^Am 2 A + O = 0.017 eV. As expected, I only 

find models predicting normal mass hierarchies below this bound. Future experiments will constrain 
many of these symmetry structures. Cosmological observations aimed at measuring the observable 
£ seem to have the best prospects. They have the potential to probe all of the inverted hierarchy 
models as well as many normal ones. Additionally, a measurement of the neutrino hierarchy could 
have a great impact on flavor symmetry models. Next generation neutrino oscillation experiments 
are expected to provide non-trivial information regarding the mass spectra. The majority are based 
on neutrino/anti-neutrino asymmetries via Earth matter effects EI ^ J _3 3 _21, 51, 52], but depend 
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strongly on large Or values. The small Or scenario is explored in [21|, 
and non-oscillation searches. 

It is important that these results and predictions are interpreted correctly. The nature of the 
fitting procedure naturally selects the lowest e value allowed by the data in order to obtain the 
smallest possible m ee . With this in mind, at least one predicted parameter should sit at the edge of 
its allowed region for each model. This does not mean that all of the explored symmetry structures 
are on the verge of exclusion. I point out that these are only predictions for the matrix that yield 
the smallest m ee value, and which is also consistent with data. If future experimental constraints 
tighten, one must simply redo the fit with the new data. Generally, larger symmetry breaking 
would be needed which would push e, and by extention m ee , higher. In terms of the histograms of 
Figure [TJ the net result would be a general movement of m ee upward. The lowest values of m ee are 
typically suppressed by the highest powers of e and will be affected the most in this transformation. 
It would lead to a narrower distribution. In this process, some models may be excluded, but it is 
not guaranteed. 

In Figure IH I show the movement of the numerical results as a function of parameter uncertainty, 
taken as a fraction of the current la deviation, as given in Table [B Here, / is defined by 

0/ = f (J current- (11.10) 

The reduced is used as input in the minimization of Eq. (|II.8|) and simulates future improvements 
in neutrino parameter measurements. This is shown in a logio scale down to / = 0.01. The black 
curves are obtained by varying the uncertainty on all of the oscillation parameters. The solid and 
dashed lines use central values equal to the best fit parameters of Table [J and equal to an alternative 
68% allowed point, respectively. The latter was chosen as a logical possibility with large deviations 
from maximal atmospheric mixing (siuOa = 0.661) and vanishing 9r {svhOr = 0.145). All other 
parameters were held at the current best fit. The solid colored curves were obtained by varying 
the parameter uncertainties individually, assuming the current best fit central values. The upper 
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FIG. 4: Plot of minimization results as a function of the fractional neutrino oscillation parameter uncertainty 
within the c = 0.5 case. The upper panel shows the variation in the smallest extracted m ee values while 
the lower panel depicts the variation in the number of allowed models. The black curves let all parameter 
uncertainties vary about the current best fit points (solid) and alternative central point (dashed) . This point 
is defined by sin0/? = 0.145 and sin 6a = 0.661 with all other parameters held at their best fit values. This 
choice is consistent with data at the 68% level. For the solid colored curves, only one uncertainty is allowed 
to vary at a time about the best fit point with the others held at their current values. 

and lower panels show the variation of the smallest derived m ee values and the total number of 
allowed models respectively with parameter uncertainty. The general trend is as expected. As the 
uncertainty is decreased, m ee is pushed larger while the number of allowed models decrease. In 
terms of individual parameter curves, sin^, and to a much lesser extent sinfl^, show the largest 
variations. All others induce very little deviation and almost sit directly on top of each other in 
the lower panel. Thus, it is clear that improved measurements of Or are essential to bounding m ee 
within this framework. 

The slope changing curve features of Figure [4] are not numerical artifacts. In the upper panel, 
they correspond to changes in the charge structure of the minimal model, each of which has its 
own m ee vs. / slope. Taking the solid black total variation curve as an example, I find that the 
six largest slope changes all correspond to optimal charge assignment changes. The exact models 
involved here are not very enlightening. As / decreases, the optimal model class jumps from C4 to 
C2 to C3. If the parameter best fit points remain at their current level, the bound on the minimum 
m ee value will triple, provided a two order of magnitude improvement in oscillation measurements. 
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Even these enhanced m ee values are beyond future experimental prospect, so the more significant 
result is related to the number of allowed models. In this scenario, the number of allowed broken 
U(l)f models is reduced by roughly a factor of three to 27. Such a small number of possibilities 
would constrain neutrino mass ultraviolet completions and help guide model builders in further 
constructions. Of course, this all depends on the exact neutrino parameter values. For example, 
the variations are not so drastic under the alternative best fit point. This is so because the majority 
of the smallest m ee structures predict large 0r and deviations from maximal a due to symmetry 
breaking effects. This point may be taken to an extreme by postulating a best fit point directly at 
the predictions of the minimal model; in which case, the lowest m ee will not vary with /, although 
the available model space might. It is not clear if the opposite extreme exists. Namely, is there a 
best fit choice that can push m ee into the next generation reach or that will narrow down the model 
space to a single, or even zero, symmetry structure? Based on current results, the former case 
seems doubtful, but the latter model reduction remains a serious possibility. This speculation may 
be verified by a scan over possible best fit parameter points, but such an analysis is not warranted 
without experimental direction in terms of improved m ee bounds and parameter measurements. 



III. CONCLUDING REMARKS 

While exact symmetries generating m ee = are excluded by current data, slightly broken 
symmetry structures are still allowed that result in m ee well beyond the reach of future (3(30u 
experiments. I systematically study these broken Abelian symmetries via a spurion analysis to 
obtain the smallest allowed m ee values. I find that, allowing reasonable 0(1) coupling constants, 
many models are excluded by data and that the smallest m ee is constrained to be larger than about 
3.9 x 10~ 6 eV at 99% confidence. The structures yielding the smallest allowed values imply general 
predictions for neutrino experiments, including large deviations from vanishing Or and maximal a, 
in addition to the normal mass spectra. Improvements in future data could increase these limits 
and perhaps, depending on the central parameter values, single out a handful of allowed models 
that may be explored in more depth. Qualitatively, the main result is that there is currently a 
small but non-zero lower bound to TppQ V that may be improved by future precision measurements. 

At face value, these results are only valid under specific circumstances. Namely, I explore models 
containing three light Majorana neutrinos subject to a broken Abelian flavor symmetry with a non- 
trivial v e transformation in the absence of fine tuning. These happen to predict small m ee bounds. 
How far can this be pushed? Abelian symmetries are better suited to this task than non- Abelian 



ones, but one may still consider discrete flavor groups. See for example [53] and references therein. 
The fact that all broken Abelian flavor models predicting m ee = are excluded helps motivate 
that no exact Abelian discrete symmetry (a discrete subgroup of U(l)f) can be allowed. This 
says nothing of non- Abelian discrete groups, which is therefore a limitation of the analysis. Next, 
relaxing the naturalness requirement, allowing the a a p parameters to take on any value, will push 
down the allowed m ee bounds, but will not admit solutions with vanishing TppQ U . The introduction 
of new degrees of freedom is more complicated, but potentially testable by other means depending 
on the masses and couplings involved. New physics above the ~ 5 TeV scale will decouple from the 
system and can be ignored for these purposes [T^]. In the intermediate range between 1 GeV and 
5 TeV, heavy particle mediation of /3(30u via contact effective operators can become important and 
influence the relationship between m ee and the effective m e JJ measured in j3j3Qv. These quantities 
should converge for small ni eej as motivated in Section [H but the detailed rate depends on the new 
physics model. Here, I only consider the mass matrix element, which will deviate from raJJ at 
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sufficiently large m ee . Hence, in the presence of TeV scale new physics, the large m ee models of 
Figure Q] may not be mapped onto Tppo„ in the standard way. Less caution is needed with the 
smallest m ee value models, which comprise the main object of this study. Moving down in mass 
scale below approximately 100 MeV, the characteristic f3f30v momentum transfer, new physics 
can mediate the decay and interfere with the light neutrinos to suppress VppQ V , regardless of the 



derived m ee value [54], |55j, |56|, l57(] . This was discussed in [55J] in terms of the eV scale type-I Seesaw 
mechanism, where the effect is particularly clear and can yield vanishing (3[3Qv rates. In this case, 
for n singlet neutrinos, the extended 3 + n x 3 + n mass matrix is relevant to the (3f3Qv system 
and the upper 3x3 diagonal block, which includes m ee , is zero by construction. Thus, the relation 
between the light neutrino mef ht defined in the last entry of Table U and r^Qv is completely spoiled. 
However, one could hope to observe these light degrees of freedom in sterile neutrino oscillation 
searches or astrophysical phenomenon. The last consideration is the presence of exact m ee = 
at some high scale where neutrino masses are generated. In the absence of a flavor symmetry to 
protect it, a non-zero m ee will be generated via two loop renormalization group effects, as motivated 
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581 ] . Details of this generation mechanism depend on new, intermediately scaled physics and has 



yet to be calculated for the e — e element of the mass matrix. Therefore, the qualitative existence of 
a lower Fppo u bound seems robust under the introduction of arbitrary new physics at scales greater 
than 100 MeV. This offers a much needed handle, or at least a conceptual proof of principle, on 
means of selecting the Dirac neutrino nature. 

This analysis is best suited to describe the theoretical possibilities of a scenario in the distant 
future where null f5f30v results push m ee bounds below the anarchy limit at 5 x 10~ 3 eV. Above 
this value, small Tpp au could be the result of mass matrix statistical fluctuations. Thus, subject 
to the above testable qualifications, it would be safe to assume that one of the flavor symmetries 
explored here is at work to suppress m ee or that neutrinos are Dirac particles. Detailed speculations 
on such broad experimental improvements are beyond the scope of this work, but it is reasonable 
that these strong f5f50u bounds would be accompanied by similar enhancements in other neutrino 
related parameters. In such a world, this analysis would guide model builders toward the complete 
and correct neutrino mass model. Currently, m ee bounds are over two orders of magnitude away 
from this situation. Still, until a positive /3/30^ signal is detected, this scenario remains a logical 
possibility that should be explored to properly understand the options open to nature. In the 
meantime, experimental bounds may be used to constrain the flavor symmetry model space. 
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